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Abstract 

The  effects  of  the  reformer  in  an  internal-reforming  molten  carbonate  fuel  cell  (IR-MCFC)  are  studied  by  mathematical  modeling. 
Temperature  distributions,  conversion  of  methane  and  compositions  of  gases  are  analyzed  through  mathematical  modeling  of  the  reformer 
and  the  cell.  In  the  reformer,  the  methane-reforming  reaction  and  the  water-gas  shift  reaction  occur  simultaneously  and  the  conversion  of 
methane  to  hydrogen,  calculated  including  the  thermodynamic  equilibrium  of  the  reaction,  reaches  99%.  Additionally,  the  endothermic- 
reforming  reaction  contributes  to  a  uniform  temperature  distribution.  The  voltage  and  the  power  of  the  IR-MCFC  are  similar  to  those  of  an 
external-reforming  molten  carbonate  fuel  cell  (ER-MCFC),  when  the  compositions  at  the  inlet  of  the  ER-MCFC  are  set  as  those  at  the  outlet 
of  the  reformer  in  IR-MCFC.  As  the  molar  ratio  of  methane  to  water-gas  decreases  at  a  fixed  total  flow  rate,  the  working  voltage  decreases. 
©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

With  regulation  of  exhaust  emissions  and  increase  in  the 
use  of  electrical  energy,  investments  in  fuel  cells  are  increas¬ 
ing.  Fuel  cells  are  of  interest  because  they  do  not  produce 
nitric  oxide  which  can  be  a  source  of  air  pollution,  have 
excellent  efficiency  in  converting  chemical  energy  into 
electrical  energy,  and  can  be  scaled-up  inexpensively. 
Molten  carbonate  fuel  cells  (MCFCs)  do  have  demerits, 
however,  in  that  they  work  at  a  high  temperature  and  the 
molten  carbonate  is  very  corrosive. 

The  temperature  distribution  in  fuel  cells  becomes  non- 
uniform  due  to  the  heat  produced  by  the  electrochemical 
reaction,  the  ion  transfer  in  an  electrolyte  plate,  and  the 
resistance  to  electron  movement  in  the  electrodes.  This 
produces  a  non-uniform  electrochemical  reaction,  thermal 
stress,  vaporization  of  electrolyte,  and  corrosion  of  materi¬ 
als.  As  a  result,  the  performance  of  fuel  cells  decreases 
rapidly.  Thus,  it  is  important  to  obtain  a  uniform  temperature 
distribution  without  locally  heated  areas  [1], 
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An  electrochemical  reaction  and  a  methane-reforming 
reaction  occur  simultaneously  in  the  internal-reforming 
MCFC.  Hence,  the  thermal  efficiency  of  the  device  increases 
and  the  system  becomes  simple,  because  the  heat  of  gen¬ 
eration  from  the  electrochemical  reaction  can  be  used 
directly  for  the  endothermic-reforming  reaction  without 
being  removed  to  the  surroundings  [2]. 

In  the  internal-reforming  MCFC,  H2  and  CO  are  formed 
by  the  water-gas  shift  reaction  (CO  +  H20  <-+  C02  +  Hi) 
and  the  methane-reforming  reaction  (CH4  +  HiO  CO  + 
3H2).  Furthermore,  the  fractions  of  reactants  and  products 
are  determined  by  the  rates  of  these  reactions,  as  well  as  by 
the  thermodynamic  equilibrium  [3 J .  Hence,  it  is  necessary  to 
consider  the  thermodynamic  equilibrium  carefully.  The  cell 
temperature,  650  °C,  is  not  sufficiently  high  to  reach  equili¬ 
brium  for  the  methane-reforming  reaction.  Therefore,  the 
fraction  of  methane  reformed  due  to  the  thermodynamic 
equilibrium  is  only  85%  [4].  Nevertheless,  the  produced 
hydrogen  and  carbon  monoxide  reach  the  equilibrium  of  the 
water-gas  shift  reaction  immediately  in  the  cell.  Hence,  the 
equilibrium  moves  in  the  direction  of  producing  hydrogen. 
So  methane  is  consumed  and  hydrogen  is  produced  by  the 
methane-reforming  reaction  and  the  water-gas  shift  reaction. 
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Nomenclature 

B  thickness  (cm) 

Cp  heat  capacity  (J/g  mol  K) 

dct  diameter  of  a  cylindrical  catalyst  (cm) 

dm  thickness  of  metal  part  in  reformer  (cm) 

E  open  circuit  voltage  (V) 

fj  material  volume  per  reformer  volume 

F  Faraday’s  constant  (96,501  C/eq.) 

(26.8  Ah/g  mol) 

h  heat  transfer  coefficient  (J/cnr  h  K) 

Hc ,  length  of  catalyst  (cm) 

I  current  density  (mA/cm2) 

k  thermal  conductivity  (J/cm  h  K) 

kr  methane-reforming  reaction  rate  constant 

(g  mol/atm  h) 

Ks  equilibrium  constant  of  the  water-gas  shift 

reaction 

L  length  of  the  fuel  cell  (cm) 

m  molar  flow  rate  (g  mol/cm  h) 

iij  mole  flux  of  fuel  gas  (g  mol/h) 

Nc  number  of  catalyst  in  Ax 

p  pressure  (bar) 

Q  heat  of  generation  per  unit  cell  area 

(J/cm2) 

R  gas  constant  (82.06  cm3  atm/g  mol  K) 

(8.314  J/g  mol  K) 

7?ohm  ohmic  cell  resistance  (Q) 

Rr  reforming  reaction  rate  (g  mol/h) 

Rs  rate  of  water-gas  shift  reaction  (g  mol/h) 

T  temperature  (K) 

V  cell  operating  potential  (V) 

W  width  of  the  fuel  cell  (cm) 

Xj  mole  fraction  of  the  zth  gas  component 

Z  effective  cell  resistance  (Q/cm2) 

Greek  symbols 

)/a  anodic  overpotential  (V) 

)/c  cathodic  overpotential  (V) 

Subscripts 
a  anode 

c  cathode 

ct  catalyst 

e  electrode-electrolyte 

i  value  at  the  inlet  of  the  gas 

m  metal 

mg  methane  gas 

re  reformer 

sg  sutTounding  gas 

up  upper  plate 

us  upper  separator 

Superscript 

0  value  at  the  entrance 


In  this  paper,  mathematical  modeling  of  a  100-cm2  fuel 
cell  is  under  taken.  The  temperature  distribution,  voltage 
distribution,  and  performance  of  an  internal-reforming  fuel 
cell  at  a  constant  current  density  are  calculated  and  com¬ 
pared  with  those  of  an  external-reforming  fuel  cell. 

2.  Theory 

A  schematic  diagram  of  the  internal-reforming  MCFC 
used  in  the  mathematical  modeling  is  shown  in  Fig.  1.  The 
reforming  section  is  composed  of  an  upper  plate,  a  catalyst 
layer,  and  an  upper  separator.  In  the  catalyst  layer,  cylind¬ 
rical  catalysts  of  radius,  r ,  and  length,  b.  are  arranged  in 
zigzag  along  the  direction  of  the  gas  flow.  In  our  mathema¬ 
tical  modeling,  however,  it  is  assumed  that  catalyst,  metal 
and  the  fuel  gas  path  are  distributed  uniformly  in  the  catalyst 
layer  with  proportions  0.22,  0.04  and  0.74,  respectively. 

Reactant  gases  flow  in  parallel  with  the  axes  of  the 
cylindrical  catalysts  with  a  velocity  v  per  unit  cross-sectional 
area.  A  Ni-Mo  catalyst  was  used.  Diffusion  into  the  pores  of 
the  cylindrical  catalysts  is  ignored.  Fuel  gas,  which  passes 
through  the  catalyst  layer  of  the  fuel  cell,  makes  a  U-turn  and 
is  used  as  the  anode  gas. 

Gases  flow  in  a  cross  type,  i.e.  anode  gas  flows  in  the 
x-direction  and  cathode  gas  in  the  y-direction.  The  paths  of 
anode  and  cathode  gases  form  a  channel.  The  channels  in  the 
gas  paths  are  neglected,  however,  and  the  paths  are  assumed 

(a)  (A)  upper  plate 


Fig.  1.  Schematic  diagram  of  (a)  internal-reforming  MCFC;  (b)  methane 
reformer.  (A)  Upper  plate,  (B)  catalyst  bed,  (C)  upper  separator,  (D)  anode 
gas  channel,  (E)  electrode-electrolyte  plate,  (F)  cathode  gas  channel  and 
(G)  lower  separator. 
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to  be  rectangular.  Additionally,  since  the  thicknesses  of  the 
separators,  the  paths  of  gases,  and  the  electrode-electrolyte 
plate  are  small  compared  to  the  lengths  and  the  widths, 
temperature  and  concentration  gradients  in  the  direction  of 
depth  are  not  considered.  The  fuel  cell  is  assumed  to  operate 
at  a  steady  state  and  the  surroundings  are  assumed  to  be 
adiabatic.  Radiation  energy  absorbed  by  gases  between  a 
separator  plate  and  an  electrode-electrolyte  plate  is  also 
neglected  [5], 

The  rate  of  the  methane-reforming  reaction  (Rr)  and  the 
rate  of  the  water-gas  shift  reaction  ( Rs)  from  an  equilibrium 
equation  can  be  obtained  from  the  following  equations, 
respectively  [6,7]: 

Rr  =  krPC  H4  (1) 

K  _  (nXco2  -fmgRs)(nXn2  +  3 Rr  - fmgRs ) 
s  _  (nXco  +  Rr  +fmSRs)(nXH2o  +  Rr  +fmgRs) 


Energy  balance  equations  for  the  methane-reforming 
section  are  as  follows.  Upper  plate  (A  in  Fig.  1): 

^sg,up  lT  T  N  ./ragging, up  , T  T  N 

D  I.  liuP  iSgJ+  „  ,  t^up  lmg) 
^up^up  ^'up'nip 

+  ^(Tup-rm)+^(rup-rct) 

^up  ^up 

+  '%%up(7-up_rus)  (6) 

^up^up 

The  boundary  conditions  are  as  follows: 

f)T 

at  x  =  0  and  x  =  L,  for  all  y  :  up  =  0  (7) 

ox 

atallx,y=0  and  v  =  W  :  =  0  (8) 

ay 


dTup  d2Tup 
dx2  dy2 


Here,  kr  (g  mol/atm  h)  is  the  rate  constant  of  the  methane- 
reforming  reaction  and  is  0.0128  at  923  K.  n  is  the  molar  flux 
of  entering  gas.  Using  Rr  and  Rs  calculated  from  Eqs.  ( 1 )  and 
(2),  the  mass  balance  equation  about  the  methane-reforming 
section  is  obtained  as  follows: 

^inVk.in  UoutXkiOUt  T  Rr  4  fmgRs  =  0  (3) 

Here,  the  subscript  k  represents  each  component  of  the  fuel 
gases,  i.e.  methane,  hydrogen,  carbon  dioxide,  steam,  and 
carbon  monoxide.  As  shown  in  Fig.  1,  methane  and  steam 
are  supplied  as  a  fuel  gas  into  the  catalyst  layer  and  produce 
carbon  monoxide  and  hydrogen.  Then,  carbon  monoxide 
and  hydrogen  are  introduced  as  an  anode  gas.  The  term  fmg  in 
Eq.  (3)  is  the  volume  fraction  of  fuel  gas  in  the  reforming 
section. 

The  same  mass-balance  equations  of  anode  gas  and 
cathode  gas  are  used  in  [5].  The  current  density,  /,  and 
the  voltage,  V,  of  cells  has  the  following  relationship  [8]: 

V  =  (FcN  —  V3n)  —  IRohm  +  Pc  ~  '7a  (4) 

Here,  the  voltages  of  the  cathode  and  the  anode,  UcN  and  VaN, 
are  obtained  from  the  Nernst  equation  and  the  overpotentials 
t]c  and  i]a  are  obtained  from  multiplication  of  current  density 
and  polar  resistance  (Zc  and  Za)  [8].  The  values  of  the  polar 
resistances  are  calculated  using  the  functions  for  the  gas 
compositions  and  the  reaction  temperature  proposed  by 
Selman  [9]  and  the  value  of  ohmic  cell  resistance,  R0 hm, 
is  calculated  using  a  function  of  temperature  [10].  By 
inserting  the  equations  for  Ecn,  Kn,  tic,  th  ar,d  Rohm,  in 
Eq.  (4),  the  following  equation  can  be  obtained  and  the  value 
of  the  voltage  can  be  obtained: 

y  =  (ya°  — yc°)--in  (x*x°2xA-iz  (5) 

v  a  cJ  2F  y  XwXaC  j 


where  h  (J/cm2  h  K)  is  a  heat  transfer  coefficient;  k 
(J/cm  h  K),  a  thermal  conductivity;  fmg  and  fm  are  volume 
fractions  of  gases  and  metal  parts  in  the  reforming  section, 
respectively.  In  the  energy  balance  for  the  upper  plate 
(Eq.  (6)),  the  terms  on  the  left-hand  side  represent  the 
conductive  heat  transfer,  and  the  terms  on  the  right-hand 
side  are  the  convective  heat  transfer  with  surrounding  gases 
and  methane  gas,  the  conductive  heat  transfer  with  metals 
and  catalysts,  and  the  radiation  heat  transfer  with  the 
separator. 

For  methane  gas  in  the  catalyst  layer  (B  in  Fig.  1): 
d(mmgCpTmg) 

—  'tup,mgv^up  -‘mg/  i  'nis,mgV^us  ^mg/ 

+  77 ie  GjCpTmg 

NcndctLcthct,mg{Tct  ) 

WAxfmg 

2NcLmBKhmmg(Tm  —  Tmg)  Qr  +  Qs 
W  Axfmg  WAxfmg 

(9) 

Boundary  condition: 

at  x  =  0,  for  all  y  :  Tmg  =  Tj  (10) 

where  Nc,  dc{  and  Lc,  are  the  number  of  catalyst  in  Ax,  the 
diameter  and  the  length  of  the  cylindrical  catalysts,  respec¬ 
tively.  <2r  and  Qs  are  the  heats  of  the  methane-reforming 
reaction  and  the  water-gas  shift  reaction. 

Metal  part  in  the  catalyst  layer  (B  in  Fig.  1): 

4NcLmBrekm(Tct  -  Tm )  |  4AcLmafmkm(Tus  -  Tm) 

dm  Bie 

.  4NcLmdmkm(Tuv  —  Tm) 

x  2NcLmBlehm^m(Tmg-Tm)  +  c  m  up - n-ll  =  0 

t>rc 


Here,  XcC  and  AaC  are  the  mole  fractions  of  CO2  in  the 
cathode  gas  and  anode  gas,  respectively.  Xw  is  the  mole 
fraction  of  H2Q  in  the  anode  gas. 
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Catalyst  part  in  the  catalyst  layer  (B  in  Fig.  1): 
^NcLmBrekcx(Tm  —  Tc  t)  4iVcLmc/m&ct(7'us  —  7A) 


2NcLmdmkct(Tup  -  Tct) 


B„ 


+  Nc7idcthmgCt(Tmg  Tc  t)  —  0 


(12) 


9  rus  |  ()  7US  /mg^mg.us  frp  rj,  ,  ,  /mg^up.us  g7.  7,  4 

v^us  -*mg/  i  ^  (-*us  — -*upj 


Upper  separator  (C  in  Fig.  1): 

92FUS  d2i 
Ox2  9v2 


Bek,,. 


Bek„e 


+  2jf(Tus-Tm)  +  2^(Tus-Tct) 


+  ^(rus-Ze)+-^(Zus-rg)  (13) 

&  s^US  ^s^ois 


where  <7m  is  the  thickness  of  the  metal  part,  as  shown  in  Fig.  1 . 

For  the  cell  section,  the  previously  used  energy  balance 
equations  are  applied  [1],  The  heat  of  electrochemical 
reaction  can  be  calculated  as  follows  [7]: 


qe  =  i 


f-TAS 
(  IF 


+  (E  —V) 


A  He  \ 
2  F+VJ 


(14) 


where  S  (J/g  mol  K)  is  the  entropy  change  due  to  the  reaction 
and  He  (J/g  mol)  the  heat  of  generation  by  the  electrochemical 
reaction  calculated  as  a  function  of  temperature,  T  (K)  [11]. 

A He  =  -57, 018  -  2.738Z  +  (0.474  x  KT3)Z2 

+  (2.637  x  KT7)Z3  (15) 


Because  both  the  methane-reforming  reaction  and  the 
water-gas  shift  reaction  occur  simultaneously  in  the  refor¬ 
mer,  the  heats  of  these  reactions  must  be  calculated.  The 
amount  of  endothermic  heat  from  the  methane-reforming 
reaction  is  obtained  as  follows: 

Gr  =  A«CH4A77r  (16) 

where  A«ch4  (g  mol/h)  is  the  number  of  moles  of  methane 
consumed  and  A Hx  (J/g  mol)  the  endothermic  heat  calcu¬ 
lated  from  the  following  relationship  [11]: 


AH,  =  190,265  +  8.314)  7.951Z- 4.354  x  4.354~3Z2 
0.097  x  105 ' 


-  3.027  x  10“6r3  - 


(17) 


The  amount  of  exothermic  heat  from  the  water-gas  shift 
reaction  is  obtained  from 


Qs  =  A«C0  AHs  (18) 

where  Anco  (g  mol/h)  is  the  number  of  moles  of  carbon 
monoxide  consumed  by  the  water-gas  shift  reaction  and  A//s 
(J/g  mol)  the  exothermic  heat  of  the  water-gas  shift  reaction 
[11],  i.e. 

A Hs  =  -9932.5  -  0.515Z+  (3.117  x  10-3)Z2 

-  (1.05  x  10"6)Z3  (19) 


Table  1 

Dimensions  of  internal-reforming  MCFC  used  in  mathematical  modeling 


Cell 

Length  (L)  (cm) 

to 

Width  (W)  (cm) 

Thickness  (cm) 

10 

Upper  separator  (5US) 

0.23 

Gas  channel  ( Bg ) 

0.2 

Anode  electrode  ( Bea ) 

0.07 

Electrolyte  plate  (Bem) 

0.1 

Cathode  electrode  (Bec) 

0.06 

Reformer 

Thickness  (cm) 

Upper  plate  (#up) 

0.23 

Catalyst  bed  (Bre) 

Catalyst 

0.12 

Diameter  (dct)  (cm) 

0.12 

Length  (Lct)  (cm) 

0.5 

Number  ( N) 

462 

Table  2 


Flow  rate  and  composition  of  fuel  gas 


Flow  rate  (mol/h) 

Composition 

Fuel  gas 

i 

ch4 

0.33 

H,0 

0.67 

The  variables  used  in  the  modeling  are  listed  in  Tables  1 
and  2.  The  area  of  the  electrode  is  100  cm2  and  all  calcula¬ 
tions  are  made  at  a  constant  current  density  of  150  mA/cnr. 
The  temperature  of  the  separator  is  923  K,  and  the  tempera¬ 
ture  of  the  entering  fuel  gas  and  that  of  the  entering  cathode 
gas  are  both  fixed  at  813  K. 

Equations  of  mass  and  energy  balances  are  calculated  by 
the  finite  difference  method.  The  operating  voltage,  tem¬ 
perature  distribution,  and  the  power  are  calculated  from  the 
modeling  of  the  fuel  cell  at  a  steady  state. 


3.  Results  and  discussion 

The  temperature  distribution,  the  cell  voltage  and  the 
performance  of  the  cell  are  predicted  by  mathematical 
modeling  of  an  internal-reforming  100-cm2  unit  cell  at  a 
constant  current  density.  The  conversion  and  composition  of 
the  gases  are  calculated.  Additionally,  changes  of  voltage 
and  power  with  current  density  are  also  determined.  The 
composition  and  temperature  distribution  are  studied  in 
order  to  evaluate  the  influence  of  the  reformer  on  the 
performance  of  the  cell.  In  addition,  the  calculated  results 
are  compared  with  those  for  an  external-reforming  unit  cell 
of  the  same  area. 

3.1.  Temperature  distribution  and  cell  voltage 

The  temperature  distribution  in  the  unit  cell  at  a  constant 
current  density  is  calculated.  The  temperature  distribution 
of  the  fuel  gas  is  shown  in  Fig.  2.  Since  the  injection 
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Fig.  2.  Three-dimensional  graph  of  temperature  distribution  of  fuel  gas  in 
reformer  section.  (B)  in  Fig.  1(b),  with  reforming  reaction  and  water-gas 
shift  reaction  at  150  mA/cm2. 


Fig.  3.  Three-dimensional  graph  of  temperature  distribution  of  anode  gas 
which  comes  from  reformer  section.  (D)  in  Fig.  1(a),  at  150  mA/cm2. 


temperature  of  the  fuel  gas  is  813  K,  the  temperature  of  the 
fuel  gas  increases  rapidly  at  the  entrance  of  the  cell  and  then 
becomes  stable  along  the  gas  flow.  Since  the  upper  plate  of 
the  reformer  is  in  contact  with  the  fuel  gas,  its  temperature 
increases  along  the  direction  of  the  fuel  gas. 

As  the  temperature  of  the  fuel  gas  increases  at  the 
entrance,  the  reforming  reaction  and  the  water-gas  shift 
reaction  are  activated  together.  The  temperature  gradient 
is  reduced  along  the  direction  of  the  gas  flow  due  to  the  heat 
of  the  endothermic  methane-reforming  reaction  and  the  heat 
of  the  exothermic  water-gas  shift  reaction.  The  heat  for  an 
endothermic  reaction  is  supplied  by  the  exothermic  reaction. 
The  temperature  distributions  become  even  because  heat 
from  the  endothermic  methane-reforming  reaction  is  more 
than  that  from  the  exothermic  water-gas  shift  reaction. 

The  temperature  distribution  of  the  anode  gas,  which  is 
injected  after  passing  through  the  reformer,  is  presented  in 


Fig.  3.  The  water-gas  shift  reaction  occurs  in  the  anode  gas. 
Since  the  cell  is  a  cross-flow  type,  the  temperature  of  the 
anode  gas  increases  in  the  direction  of  the  cathode  gas  flow.  It 
decreases  in  the  direction  of  the  anode  gas,  however,  because 
of  the  low  temperature  of  the  fuel  gas  in  the  reformer. 

The  temperature  distribution  of  the  electrode-electrolyte 
plate  is  shown  in  Fig.  4.  It  is  influenced  by  the  anode  and  the 
cathode  gases  and  varies  in  the  direction  of  the  cathode  gas. 
The  temperature  distribution  of  the  electrode-electrolyte 
plate  changes  under  the  influence  of  the  cathode  gas.  The 
temperature  distribution  is  quite  even  in  the  direction  of 
the  anode  gas. 

The  temperature  distribution  of  the  electrode-electrolyte 
plate  is  compared  in  Fig.  5  for  MCFCs  with  and  without 
internal  reformers.  There  is  little  difference  in  the  distribu¬ 
tion  in  the  direction  of  the  cathode  gas  flow,  as  shown  in 
Fig.  5(b).  As  shown  in  Fig.  5(a),  however,  the  temperature 


Fig.  4.  Three-dimensional  graph  of  temperature  distribution  of  electrode-electrolyte  plate,  (E)  in  Fig.  1(a),  at  150  mA/cm2. 
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Dimensionless  distance 
(a)  Anode  gas  *4 


Dimensionless  distance 
(b)  cathode  gas 

Fig.  5.  Temperature  distribution  of  electrode-electrolyte  plate  along 
direction  of  (a)  anode  gas  and  (b)  cathode  gas. 


distribution  in  the  direction  of  the  anode  gas  flow  become 
uneven  due  to  the  internal  reformer.  This  is  because,  the  heat 
from  the  electrochemical  reaction  is  consumed  by  the 
methane-reforming  reaction  in  the  reformer. 

3.2.  Conversion  of  methane 

The  main  component  of  the  fuel  gas  is  methane  which 
changes  into  hydrogen  in  the  reforming  section.  Here,  there 
are  two  reactions  of  fuel  gas,  namely,  the  endothermic  methane 
reaction  and  the  exothermic  water-gas  shift  reaction.  These 
two  reactions  affect  cell  performance  and  temperature  dis¬ 
tribution.  Additionally,  hydrogen  produced  by  the  methane¬ 
reforming  reaction  is  used  as  the  anode  gas.  Thus,  the  effects 
of  the  methane-reforming  reaction  are  analyzed. 

The  conversion  of  methane  and  composition  of  each 
component  in  the  direction  of  the  gas  flow  are  shown  in 
Figs.  6  and  7,  respectively.  The  methane-reforming  reaction 
depends  on  the  concentration  and  temperature.  Hence,  due 
to  the  rapid  rise  of  temperature  at  the  entrance  (Fig.  5),  the 
methane  conversion  becomes  0.4  rapidly  at  the  entrance  and 


Fig.  6.  Changes  of  conversion  of  fuel  gas  in  reformer. 


then  rises  gradually  to  1  in  the  reforming  reaction  and  the 
shift  reaction  while  the  fuel  gas  passes  through  the  reformer, 
as  shown  in  Fig.  6. 

At  the  entrance  of  the  reformer,  the  concentration  of 
methane  is  very  high  and  the  reaction  in  the  catalyst  is  very 
active.  As  the  amount  of  methane  decreases  in  the  direction 
of  the  fuel  gas,  it  decreases  gradually.  Hence,  as  shown  in 
Fig.  7,  the  amounts  of  H20  and  methane  decrease  and  the 
amounts  of  hydrogen  and  carbon  monoxide  increase.  The 
amount  of  hydrogen  increases  by  a  factor  of  three  compared 
with  the  decrease  in  the  amount  of  methane.  Hydrogen 
produced  in  the  reformer  is  used  as  a  fuel  gas.  The  total 
flow  rate  of  the  fuel  gas  in  the  reformer  increases  as  the 
reactions  proceed.  The  conversion  of  methane  in  the  refor¬ 
mer  is  more  than  99%. 

The  variation  in  average  cell  voltage  with  current  density 
at  different  ratios  of  fuel  gas,  CH4:H20,  is  shown  in  Fig.  8. 
The  current  density  is  held  constant.  The  voltage  and  power 
at  105  mA/crn2  are  found  to  be  1.05  V  and  15.75  W,  respec¬ 
tively.  The  compositions  of  the  entering  gas  change,  but  the 


Fig.  7.  Changes  in  mole  fraction  of  each  component  in  fuel  gas  at 
150  mA/cm2. 
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Fig.  8.  Changes  of  cell  voltage  with  current  density  at  different  fuel  gas 
ratios,  CFL^F^O. 


Current  density  (mA/cmz) 

Fig.  10.  Changes  of  average  cell  voltage  with  current  density  with  and 
without  reformer. 


entering  total  flow  rate  is  fixed.  The  best  performance  of  the 
cell  is  obtained  at  a  ratio  of  CrLphTO  =  1:2.  A  decrease  in 
the  ratio  at  a  fixed  total  flow  rate  means  the  amount  of 
CH4  decreases.  Hence,  the  amount  of  hydrogen  from  the 
methane-reforming  reaction  decreases.  As  the  fraction  of 
hydrogen  in  the  fuel  gas  entering  into  the  anode  increases, 
the  working  voltage  tends  to  increase.  The  working  voltage 
is  small  at  CH4:H20  =1:1.  This  is  because,  there  is  remain¬ 
ing  methane  from  the  reformer,  the  produced  hydrogen  is 
small  and  the  temperature  becomes  low.  As  a  result,  the  rate 
of  the  electrochemical  reaction  decreases,  the  resistance  to 
the  movement  of  ions  increases,  and  the  electrode  resistance 
of  the  cell  increases.  Additionally,  the  power  of  the  cell 
increases  with  increase  of  hydrogen  fraction,  since  the 
power  is  proportional  to  the  amount  of  consumed  hydrogen. 
The  calculated  voltages  at  a  constant  current  density  of 
150mA/cm2  are  given  in  Fig.  9.  There  is  a  difference  of 
0.1  V  in  the  whole  cell. 


Dimensionless  distance 
anode  gas 

Fig.  9.  Contour  graph  of  cell  voltage  distribution  at  constant  current 
density  of  150  mA/cnr. 


3.3.  Effect  of  current  density  on  cell  performances 

The  working  voltage  and  power  in  internal-  and  external- 
reforming  MCFCs  at  each  current  density  have  been 
calculated.  Changes  in  the  working  voltage  with  current 
density  are  compared  in  Fig.  10.  In  this  comparison, 
conditions  such  as  the  flow  rate  and  composition  of  the 
entering  anode  gas  in  the  external-reforming  fuel  cell  are 
matched  with  those  of  the  exiting  gas  from  the  reformer  in 
the  internal-reforming  fuel  cell.  As  shown  in  Fig.  10,  the 
voltage  depends  almost  linearly  on  the  average  current 
density.  This  is  because  the  voltage  decreases  due  to  the 
ohmic  resistance  of  the  electrolyte.  Voltages  and  current 
densities  in  the  internal-reforming  fuel  cell  are  slightly 
higher  than  those  in  the  external-reforming  fuel  cell. 
The  differences  in  temperature  distribution  are  also  very 
small. 


4.  Conclusions 

The  effects  of  the  reformer  in  the  internal-reforming 
MCFC  have  been  studied  by  mathematical  modeling.  The 
temperature  distribution,  the  conversion  of  methane  and  the 
composition  of  gases  are  analyzed  through  mathematical 
modeling  of  the  reformer  and  the  cell. 

In  the  reformer,  the  methane-reforming  reaction  and  the 
water-gas  shift  reaction  occur  simultaneously  and  the  endo¬ 
thermic-reforming  reaction  contributes  to  a  uniform  tem¬ 
perature  distribution.  The  flow  rate  increases  as  the  reactions 
proceed.  Most  of  the  methane  converts  to  hydrogen  under 
the  conditions  studied. 

The  reactants  change  rapidly  at  the  entrance  due  to  a 
temperature  rise.  Hence,  the  rate  of  the  methane-reforming 
reaction  increases.  The  equilibrium  constant  of  the  water- 
gas  shift  reaction  becomes  small  at  a  high  temperature.  As  a 
result,  the  rate  of  the  methane-reaction  increases  quite 
slowly  due  to  the  water-gas  shift  reaction. 
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As  the  current  density  increases,  the  cell  voltage  decreases 
and  the  power  increases.  On  changing  the  ratio  of  the 
fraction  of  methane  to  that  of  water-gas,  the  largest  value 
of  the  working  voltage  is  obtained  when  the  ratio  is  1:2. 

When  the  conditions  of  the  entering  anode  gas  in  the 
ER-fuel  cell  are  matched  with  those  of  the  exiting  gas  from 
the  reformer  in  the  internal-reformimg  fuel  cell,  voltages 
and  current  densities  in  the  internal-reforming  fuel  cell 
are  slightly  higher  than  those  in  the  external-reforming  fuel 
cell. 
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